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Abstract 



I We report a measurement of the rate for B — > 77 A s transitions where the 77 meson has center-of- 
■ mass momentum in the range 2.0 to 2.7 GeV/c and X s represents a system comprising a kaon and up 
to four pions. Our study is based on 22.2 million BB pairs collected at the T(4S) with the BABAR 
detector at the Stanford Linear Accelerator Center. We find B(B —> rj'X s ) = (6.8~^®'l(stat) ± 
1.0(syst)*Q'®(bkg)) x 10~ 4 assuming that the signal is due to b — > sg* transitions. 
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1 Introduction 



The study of B decay modes involving gluonic penguin transitions b — > sg* is important both for 
obtaining a better understanding of the mechanisms contributing to B decays and as a sensitive 
place to search for direct CP violation effects. Exclusive modes such as B + — > rj'K + or B° — ► K + it~ 
are expected to be dominated by penguin amplitudes .T] and may show direct CP-violating charge 
asymmetries j2]- The study of a collection of decay modes B — * rj'X s , where X s denotes a set 
of decay particles containing an s quark, is another attractive method for obtaining inclusive 
information about b — > sg* transitions. 

In this paper, we report an application of this last approach to obtain a semi-inclusive mea- 
surement of the rate for B — > rj X s for center-of-mass momentum p*, of the rj ranging from 2.0 
to 2.7GeV/c. In this high momentum interval, rj production from b — > c — > rj cascades, such 
as B D S X with D s r/X, B D + X with D + B £>°X with D° r/X, 

B — > yl c X with yl c — > r/'X, is suppressed, although it is important to note that other B de- 
cay processes, such as b — > u decays (f? — ► 7/71", r/'p, rj'ai) and internal spectator 6 — > c decays 
(5° — > r]'D (*') may still contribute. Figure ^ shows the momentum distribution for the rel- 
evant processes. The various contributions can be distinguished by kaon identification and by 
examining the mass spectrum of the system recoiling against the rj {M{X S )). This is possi- 
ble if the recoiling hadronic system is reconstructed with a semi-inclusive set of possible decay 
modes. CLEO used such a semi-inclusive reconstruction technique to obtain the branching fraction 
B(B -> rfX s ) = (6.2 ± 1.6(stat) ± 1.3(syst)±°{%(bkg)) x 10" 4 for 2.0 < p*, < 2.7 GeV/c ;3j. This 
rate is large in comparison with available predictions [3]. Our analysis uses a similar technique, 
although with improved kaon identification and different Monte Carlo models. 

2 The BABAR detector and dataset 

The data used in this analysis were collected with the BABAR detector [5] at the PEP-II storage 
ring 6 located at the Stanford Linear Accelerator Center. The study presented here is based on 
an integrated luminosity of 20.2 fb _1 corresponding to a sample of 22.2 million BB pairs at the 
T(AS) resonance (on-resonance) and 2.6 fb -1 collected with center of mass energy 40MeV below 
this resonance (off-resonance). 

The asymmetric beam configuration in the laboratory frame provides a boost to the T(4S) 
increasing the momentum range of the i?-meson decay products up to 4.3 GeV/ c. Charged particles 
are detected and their momenta are measured by a combination of a silicon vertex tracker consist- 
ing of five double-sided layers and a 40-layer drift chamber, both operating in a 1.5-T solenoidal 
magnetic field. Photons and electrons are detected by a CsI(Tl) electromagnetic calorimeter, which 
provides excellent angular and energy resolution with high efficiency for energies above 20MeV. 

Charged particle identification is provided by measurements of the average energy loss dE/dx 
in the tracking devices and the Cherenkov angle in the detector for internally reflected Cherenkov 
light (DIRC). The Cherenkov angle determination provides K-ir separation of better than 4<r below 
3GeV/c and 2.5a for the highest momenta. 

3 Event selection 

In order to select BB events, we use the following requirements: 
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Figure 1: Distribution of p*, from Monte Carlo simulation of various b — > sqq processes, such as 
B — ► r/'K, rj'K* (horizontally hatched histogram), rj from a mixture of BB sources dominated 
by b — > c cascade contributions (diagonally hatched histogram), r/'D (open histogram) and r/'D* 
(dashed histogram) samples. 



• There must be at least four charged tracks per event, in order to suppress low multiplicity 
events such as e + e~ — > e + e~, t + t~. 

• i?2, the ratio of the second to the zeroth Fox- Wolfram moment [7j, must be less than 0.5. The 
distribution of this variable, which varies from to 1, is peaked at low values for spherical 
BB events while it is broader and shifted toward intermediate values for e + e~ — > qq events 
(q = u,d,s,c). 

• The total energy of charged and neutral particles is required to be at least 5 GeV and below 
15 GeV. Events having a significant missing energy due to neutrinos, such as t + t~ events 
where one r decay to a higher multiplicity hadronic mode, such as r — ► 7nriri , T , are highly 
suppressed by this cut. 

The efficiency of this selection is (98 ± 1)% for BB Monte Carlo events. 

4 Semi- inclusive analysis 

We form combinations of a charged kaon or a Kg, an rj' and as many as four pions of which at 
most one is a 7r°. Sixteen decay modes and their charge conjugates are considered: 
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The rj' is reconstructed in the rjir + ir~' channel, from ij — > 77 candidates only. We require a 
minimum energy of 50 MeV for photons from 77 — > 77. Photons are rejected if they are consistent 
with originating from a 7r° having an energy above 200 MeV. Candidates with invariant mass 
within 3a of the 77 mass are kinematically fitted to the nominal 77 mass and then combined with 
two charged tracks to form an 7/ candidate. 

To identify the s quark in the X s system, we select events either with a track consistent with 
a charged kaon or with a reconstructed K® in the 7r + 7r~ channel. The charged kaon selection 
has been optimized to reduce background from B — > r/V, ifp, and rj a\ decays without losing too 
much efficiency for the signal. For the K®, we consider the two-dimensional angle a between the 
momentum of the K® candidate and the flight direction, which is peaked near zero for a true K®. 
We require a < 0.05 rad. Selected K® candidates with invariant mass within 3a of the K® mass 
are kinematically fitted to the nominal K® mass. 

4.1 B candidate selection 

Candidates for B — > r]'X s are reconstructed in the 16 modes listed above. They are required to be 
consistent with a B decay based on the energy substituted mass 0, 



where y/s is the total e + e~ center-of-mass energy. The initial-state four- momentum (E^, pj) derived 
from the beam kinematics and the four-momentum (Eb,Pb) of the reconstructed B candidate are 
all defined in the laboratory. The calculation of Mes only involves the three-momenta of the decay 
products, and is therefore independent of the masses assigned to them. An additional variable, the 
cosine of the angle in the center of mass frame between the thrust axis of the B candidate and the 
thrust axis of the remainder of the event, cos9t, is used to remove continuum background. The 
selection criteria applied are the following: 

• Mes > 5.265 GeV/c 2 and \AE\ < 0.1 GeV, i.e., consistent with the nominal B mass and 
known production energy; and 

• I cos 9t\ < 0.8 to reduces the large continuum background, which is concentrated near cos 9t = 
±1 while the expected signal is uniformly distributed in this variable. 

For each event, we select only a single candidate in a given decay mode. The selected candidate is 
the one with the smallest x 2 defined as x 2 = {Mes ~ M B (PDG)) 2 /cr(M|) + AE 2 /<r(AE 2 ), where 
the widths a(Mg) and a(AE 2 ) are obtained from Monte Carlo simulation. 




and the energy difference 



AE = (EiE B — pi ■ ps — s/2)/sfs 



10 



4.2 Determination of the X s mass spectrum 



To explore the raw X s mass distribution, we first select the B candidates for which the mass of the rf 
daughter is within 3 sigma of the known value, and subtract the off-resonance contribution, rescaled 
by the luminosity ratio, from the on-resonance distribution. The resultant mass distributions for 
all B modes and separately for the B° modes are shown in Fig. [2j Both distributions can be seen 
to peak above 2GeV/c 2 . We can also examine the X s mass spectrum for a possible signal for the 
internal spectator charmed decays (X s = D°(*> with D*° — ► D°ir°, D° — ► Kir, Kit-kit, Ktttt ). In 
particular, for the B° modes, there is no evidence for a narrow D° signal near 1870MeV/c 2 (see 
Fig. 01 for predicted distribution), although statistics are low. 




Figure 2: Continuum-subtracted M(X S ) spectrum for all B modes (left) and B° modes alone (right) 



To obtain the decay X s spectrum we first fit the rf mass distribution in bins of X s mass. 
For masses above 2.32GeV/c 2 , corresponding to the kinematic limit p*, < 2GeV/c, the yield is 
dominated by the b — ► c — > rf contribution. The differential branching fraction for the region 
M{X S ) < 2.5GeV/c 2 , where we expect b — > sg* to be dominant, is shown in Fig. 21 The signal 
tends to peak towards higher mass values, and remains substantial between 2 and 2.5GeV/c 2 . The 
experimental resolution for M(X S ) has been estimated with the Monte Carlo simulation to be 
80-90 MeV/c 2 around 1.5GeV/c 2 , rising to 170-180 MeV/c 2 for 2GeV/c 2 and above. 

Looking specifically at the two-body decay modes alone (X s = K^, K®), we find no significant 
signal for rj'Kg but observe 36.2 ±6.6 events for rfK^ (see Fig. EJ). This corresponds to a branching 
fraction of (5.6 ± l.O(stat)) x 10~ 5 , in good agreement with our recent exclusive measurement of 
(7.0 ± 0.8(stat) ± O.b(syst)) x 10~ 5 l^j, thereby confirming the consistency of the semi-inclusive 
method. 

4.3 Extraction of rf signal 

The B — ► rfX s yields are determined by fitting the observed rf signals, after all selection criteria 
have been applied in the range 2.0 < p*, < 2.7GeV/c. We fit separately the K ± modes and the 



11 




0.5 1 1.5 2 2.5 

M(Xs) GeV/c 2 



Figure 3: M{X S ) spectrum predicted from simulation of B — > rj'D decays 
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Figure 4: Distribution of dB(B — > 7]'X s )/dM(X s ) (statistical errors only) as a function of M(X S ) 
for all B modes. 



Kg modes with results as shown in Fig. H3 We find 188.8 ± 21.5 events for the modes and 
57.1 ± 14.7 events for the K® modes in on-resonance data. For off-resonance data, we find O.O^q'q 
for modes and O.O^q'q for K® modes. 
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Figure 5: rf mass for X s = K^. 



The efficiencies have been computed with two main Monte Carlo simulations to study model 
dependence. One involves a mixture of resonant modes, rj'K, rj'Ki, rj'K^, Tj'K^, and rf K\, and 
the other uses an X s pseudo particle decaying to sqg (q = u,d). In the latter, we use a M(X S ) 
distribution derived from the rf QCD anomaly theoretical prediction [HJ |2] . 

For a given K mode (K^ or K®), the efficiency is computed as: 



where iVgt and iVg en are, respectively, the numbers of fit and generated events. The efficiencies are 
estimated (statistical errors only) to be (7.6 ±0.4)% for K modes and (3.3 ±0.2)% for modes, 
including the Kg — > tt + tt~ branching fraction. 

4.4 Semi-inclusive branching fraction and interpretation 

The semi-inclusive rate for the high momentum region, 2.0 < p*, < 2.7GeV/c, is computed by 
performing a weighted average of the results obtained in the modes and the K® modes. The 
detection efficiencies given in Section POI are corrected to account for the rf and rj branching fractions 
to the channel we observe (17.4%), and the K° — > K® projection (1/2). The final state X s includes 
both K + - and X°-tagged decays, and both charged and neutral B mesons contribute to the observed 
yields. The total number of produced B mesons of both charges is 45.5 million. Computing the 
weighted average of the K + - and -fC°-tagged decays, assuming that their branching fractions are 
equal, we obtain B(B rj'X s ) = (6.8til(stat) ± l.0(syst)±Q%(bkg)) x 10" 4 . Our measurement 
confirms, with higher precision, the CLEO result, (6.2 ± 1.6(stat) ± 1.3(syst)ti° 5 {bkg)) x 10" 4 . 

Sources of systematic error included in this result are summarized in Table H The largest 
uncertainty arises from our ability to model the X s system. This contribution is estimated by 
determining the efficiency with different Monte Carlo generators, as described in Section [4.31 The 
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Figure 6: Semi-inclusive rf mass spectra for (upper left: on-resonance, lower left: off-resonance) 
and (upper right: on-resonance, lower right: off-resonance) modes, in the momentum range 
2.0 <p*, < 2.7GeV/c 



efficiency from the X s pseudo-particle model tends to be lower than that obtained with the mixture 
of resonant decays, the difference is smaller when we increase the amount of heavier resonances (-K3, 
if |) in the model. Other systematic uncertainties are contributed by the event shape requirement, 
BB counting, tracking and photon detection efficiencies, and kaon identification efficiencies. 

We have also considered a possible contribution from internal spectator decays 5° — > rfD ^ 
as an additional source of systematic error. Early branching fraction predictions for these modes 
were in the range (1.5 — 6.0) x 10 -5 jllj while more recent predictions |12j find (3 — 5) x 10 -5 , 
assuming that only the quark content of the rj' contributes to the tree diagram. The efficiency 
computed for these modes is 2%, which implies a contribution of (1.4-4.6)% to the observed yield. 
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Table 1: Systematic errors for B(B — > rj'X s ) 



K ± modes 



Kg modes 



Source 



Model dependency 
Event shape 
B counting 
Tracking 
Photon detection 

B(rj' — ► 7/7T7T) x £>(r? 

i^ 1 * 1 identification 



77) 



±10% 
±1% 
±1.6% 
±5% 
±6.8% 
±3.4% 
±2.6% 



±10% 
±1% 

±1.6% 
±7% 

±7.2% 

±3.4% 



To conservatively account for the theoretical uncertainty, we use a branching fraction of 9.0 x 10 
in determining the systematic error contribution. 

Several models have been proposed to explain rf production at high p*, in B decays: 

• Introduction of b — > sqq operators with constructive interference between uu, dd and ss 
components of the rj (3]; however the expected branching fraction from such processes is 
much lower than we observe. 

• A b —* ccs enhancement through a possible cc content in the rj wave function |13U14j : however, 
more recent calculations |15| do not support such a wave function component. 

Both of these mechanisms predict an X s mass spectrum peaking near 1.5GeV/c 2 , lower than sug- 
gested by the data. The branching fraction and mass spectrum in data do appear to be consistent 
with a model incorporating an rf QCD anomaly jHJ ^ that couples the rf to gluons, b — > sg*, 
9 -> 9V ■ 

5 Summary 

We have measured a semi-inclusive rj branching fraction using fully reconstructed final states 
consisting of an 77' and a system comprising one kaon and up to four pions to reduce the background 
from other B decays. We find B(B -> rfX a ) = (6.8±$(atot) ± l.O(syst)t^) x 10~ 4 for 2.0 < p*, < 
2.7GeV/c. This measurement confirms the surprisingly large branching fraction found by CLEO, 
but with much improved precision. This and the shape of the X s mass distribution are consistent 
with the 7/ QCD anomaly model. 
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